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Abstract 

o . 

A model is introduced for two reduced BCS systems which are 
coupled through the transfer of Cooper pairs between the systems. 
The model may thus be used in the analysis of the Josephson effect 
l/-) ■ arising from pair tunneling between two strongly coupled small metal- 

lic grains. At a particular coupling strength the model is integrable 
and explicit results are derived for the energy spectrum, conserved 
operators, integrals of motion, and wave function scalar products. It 
is also shown that form factors can be obtained for the calculation of 
correlation functions. Further, a connection with perturbed conformal 
& • field theory is made. 
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1. Introduction. 



The study of superconductivity in metallic nanoparticles is of fundamen- 
tal scientific interest and is the basis of several proposed solid state quantum 
computers The experimental work of Ralph, Black and Tinkham || 
has stimulated a wealth of theoretical activity concerning the reduced BCS 
Hamiltonian describing electron pairing correlations in systems with discrete 
energy levels and a fixed number of electrons. For sufficiently small systems 
the energy level spacing becomes comparable to the energy gap predicted 
by the BCS mean-field approximation making mean-field theory unreliable 
because of the large superconducting fluctuations . 

These considerations are also relevant to proposals to observe BCS su- 
perconductivity in gases of fermionic atoms such as spin-polarised 6 Li M. 
Quantum degeneracy of 6 Li at temperatures of about 240 nK has recently 
been observed in an atom trap with frequencies, u ~ 1 kHz [J5J, corresponding 
to an energy level spacing of the order of 10 -12 eV. The estimated BCS tran- 
sition temperature is of the order of 20 nK ||, corresponding to an energy 
gap of the order of 4 x 10~ 12 eV, comparable to the energy level spacing. 

The failure of approximation schemes increases the importance of the fact 
that the reduced BCS Hamiltonian is exactly solvable ||. Furthermore, it 
is integrable J7|, and has a rich mathematical structure [§, ||, |i~0f . The in- 
tegrability of the model has recently been clarified by the Quantum Inverse 
Scattering Method (QISM) and algebraic Bethe ansatz || [10|1 . This has the 
advantage that it makes accessible the computation of form factors and cor- 
relation functions [§, |9| . This realisation makes it clear that these techniques 
can be used to study more general classes of pairing models. For example, 
the coupled pairing model which accomodates proton-neutron pairing in nu- 
clear systems studied by Richardson |y]] can be formulated via the QISM 



utilizing the so(5) symmetry algebra pi. Other examples are discussed in 
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For a single metallic grain, due to the fixed and finite number of electrons 
there is no breaking of gauge symmetry and no long-range superconducting 
order. Consequently, an important question concerns the existence and na- 
ture of the Josephson effect [15] between two coupled metallic grains. If the 
number of pairs of electrons N in each grain is fixed, then for very weak 



intergrain coupling there is no Josephson effect [16|. On the other hand, if 
one constructs a qubit consisting of a coherent superposition of N and N + 1 
pairs in an individual grain, Josephson tunneling occurs and can be used to 
couple qubits together fl], [16| . 

In this article we introduce a model describing the tunneling of Cooper 
pairs between two metallic grains described by reduced BCS models. We 
show that for a particular value of the intergrain coupling the model is inte- 
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grable and has an underlying so(3) symmetry which we exploit to solve the 
model by the algebraic Bethe ansatz. First we consider the physical basis of 
the model, stressing that it requires a large intergrain coupling. 



2. Model Hamiltonian. 



The physical properties of each metallic grain are described by the reduced 
BCS Hamiltonian g| 

c c 
H B cs(q) = e i n M) - 9 Yl h l^) h M) 

with q = 1, 2. Above, j = 1, labels a shell of doubly degenerate single 
particle energy levels with energies ej and rij is the fermion number operator 
for level j. The operators bj = Cj^Cj.i and &j are the annihilation and creation 
operators for the Cooper pairs and satisfy the hard-core boson relations 

(b]Y = 0, [&„ 6+] = S jk (l - 26J&,-) 

[b j ,b k ] = [b],bl\=0, for MJ. 

One of the prominent features of this Hamiltonian is the blocking effect. For 
any unpaired electron at level j the action of the pairing interaction is zero 
since only paired electrons are scattered. This means that the Hilbert space 
can be decoupled into orthogonal subspaces of paired and unpaired electrons 
in which the action of the Hamiltonian on the subspace for the unpaired 
electrons is automatically diagonal in the natural basis. Consequently, we 
need only focus on the subspace of paired electrons when solving for the 
spectrum of the model. 

In order to firmly establish a microscopic basis for Josephson's proposal 
|TT| , |T8f of Cooper pair tunneling through an insulating barrier between two 
bulk superconductors, Anderson 16(1 , and Ambegaokar and Baratoff (19 



studied the Hamiltonian 

H = HbcsO) + Hbcs(2) - ]Tt ifc (4(1)0^(2) + C J )(7 (2)c fe , ff (l)) . (1) 

They treated the effect of the tunneling to second-order in tjk and the BCS 
Hamiltonians were treated at the mean-field level in the grand-canonical 
ensemble. If we define A exp(i0(g)) = g £\ (cj j(<?)cj ^(q)) then the Josephson 
current is proportional to t 2 Asin(0(l) — 0(2)) where t 2 = Ylijk \tjk\ 2 - One is 
then justified in writing the effective Hamiltonian 

c 

H = H BCS {l)+H BCS {2)-ejY{b){l)b k {2)+b){2)b k {l)) (2) 

j,k 
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where Ej is the Josephson coupling energy. 

If one repeats the calculation with N(l) and N(2) fixed and finite then to 
second-order in the tunneling matrix element there is no Josephson coupling 
because (ct ^{q)c^(q)) = for both q = 1,2. However, if the intergrain 
coupling is sufficiently large that we go beyond second-order in t then only 
N(l)+N(2) is a good quantum number, and not N(l) and N(2) individually. 
It is then possible that (ct^(g)ct , (q)) ^ and so it is appropriate to study a 
Hamiltonian of the form (Q). 

The Hamiltonian (§) is not integrable for arbitrary values of the cou- 
pling parameters g and Sj. We demonstrate below that for the particular 
choice Ej = g the model does become integrable and an exact solution can be 
obtained. In typical Josephson junctions involving spatially separated super- 
conductors that are separated by an insulating barrier, ej <C g- However, this 
does not rule out fabrication of a junction consisting of two metallic grains 
separated by a quantum point contact with a transmission coefficient close 
to one. Because of the tunability of the pairing interaction in trapped atom 
systems, as noted in |EJ, the case of (relatively) strong Josephson coupling 
between such systems may be feasible. These systems also admit the pos- 
sibility of an "internal" Josephson junction which does not involve spatially 
separated condensates EO . 



3. Integrable coupling model. 



As in the case of a single BCS system, the model (0) exhibits the blocking 
effect for unpaired electrons. However, the model with ej = g displays an 
additional blocking effect which we now elucidate. We begin by writing the 
Hamiltonian in an equivalent form 

c c 

H = J2 W ~ 9 £(*}(!) + 6}(2))(6*(1) + h(2)) (3) 

where rij = n,j(l) + nj(2). The pairing interaction cannot scatter any anti- 
symmetric state of the form 




since an elementary calculation shows that 

(6 fc (l) + 6 fc (2))|* i ) = 

for all k and j. A consequence of the blocking effect for antisymmetric pairs 
is that the Hilbert space now decomposes into subspaces where the pairing 
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interaction is non-trivial on the local states spanned by 

|1> = &W(2)|0>, 
|2) = -i= (&t(i) +6 t (2) ) | > , 

|3> = |0) 

for each level j. We will denote this space Vj. Hence the problem of diag- 
onalizing the Hamiltonian is reduced to diagonalizing the Hamiltonian over 
the product of these subspaces. This will be achieved below using the alge- 
braic Bethe ansatz and the spin 1 representation of the so(3) Lie algebra. 
The procedure is analogous to that employed in Ref. || [H| for the standard 
reduced BCS model involving su(2). 

The hard-core bosons provide a spin 1 representation of so(3) through 

S+= (6 t (l) + 6 t (2)), S- = (6(1) + 6(2)), S* = N - 1 (4) 

where N = N(l) + N(2) counts the number of Cooper pairs. Associated with 
the spin 1 representation there is a solution of the Yang-Baxter equation 

Ri 2 (u - v)L 1 (u)L 2 (v) = L 2 {v)Li(u)Ri 2 {u - v) 

with 

2 

r{u) = + 1 C®C 

m,n=l 

L(u) = I®I+ 1 l(e\®S z -el®S z + e 1 2 ®S- + e 2 1 ®S + ) 

where {e™} are 2x2 matrices with 1 in the (m, n) entry and zeroes eleswhere. 
Above, / is the identity operator and rj is a scaling parameter for the rapidity 
variable u which plays an important role in the subsequent analysis. With 
this solution we construct the transfer matrix 

t(u) = tr (G L 0C (u - e c )...L 01 (u - e x )) (5) 

which is an operator on the subspace V\ ® V 2 ® ---Vc on which the pairing 
interaction acts non-trivially. Above, we have introduced an auxiliary space 
Vo and tro denotes the trace taken over this space. Moreover G = exp(ar)a) 
with cr = diag(l, —1). A consequence of the Yang-Baxter equation is that 
[t(u), t(v)] = for all values of the parameters u and v. Defining 

Ta = lim - — ^-t(u) 

for j = 1, 2, C, we may write in the quasi- classical limit Tj = Tj + 0(77) and 
it follows that 

fa, T k \ = 0, Vj, k. 
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Explicitly, these operators read 




with 6 = S + <g> S + S <g> >S + + 2S Z <8> S z . It is an algebraic exercise to show 
that the Hamiltonian (^j can be expressed as 

c c c c 

3 j,k j j 

with g = —1/at and the elements of so(3) represented through (Q). This shows 
that the Hamiltonian is integrable since it is clear that [H, Tj] = 0, Vj. 



4. Exact solution. 



A standard calculation (e.g., see pip gives that the eigenvalues of the 
transfer matrix (|5]) take the form 

c M 
A(u) = exp {ari)\ I M J - 

k U ~ €k j U ~ W i 

C M 

+ exp(-a?7) I I I I J - 

V u ~ e k y u- Wj 

where the parameters Wj are required to satisfy the Bethe ansatz equations 



c M 



exp 



Wl~e k + T] j-jWl-Wj+T) 



, Wl-e k -T) Wi - Wj - T] 

K J 



The eigenvalues of the conserved operators (||) are obtained through the 
leading terms of the expansion of the transfer matrix eigenvalues in the pa- 
rameter 7]. This yields the result 

C. M 

such that the parameters Vj satisfy the coupled algebraic equations 

— + E^- = E — ( 8 ) 



which are analogous to Richardson's equations || for the reduced BCS model. 
Through (^), the energy eigenvalues are found to be 



c M 
E = AY,e J -2Y,v 3 + 2g{M-C). 

3 j 

We introduce a set of states of the form 

M 



\w, M) = \w h ...w M ) = Yl B(wj) |$) (9) 

j 

where |$) is the completely filled state with 2C Cooper pairs and 

c 

B(u) = Y, 



1 Si 



U — 6j 

1 J 



Note that these operators are mutually commuting for different values of u. 
In general we find that 

Tj \w, m) = \j \w, m)-J2 T^rtr 1^' M ^ ' ( 10 ) 

where 



Ej - W a 
a J 



2 " 2 A 2 



and Aj is understood to be a function of the parameters Wj. In ([II]) we 
defined \w, M) a by 



C S7 



\w 



^— ' 111 f: ■ 



W a - tj 



Imposing f a = 0, one immediately sees that \w, M) becomes the eigenvector 
of the integrals of motion Tj with Xj as the eigenvalue. The constraint f a = 
is then equivalent to ([8]). For a given solution of (||), the corresponding 
eigenstate has 2C — M Cooper pairs. 



5. Connection to conformal field theory. 



Introduce the function x(t>, e) = x{ v i> " " • > v n, £i, • •• , Cn) obeying the fol- 
lowing differential relations 

dx(v,e) . dx(v,e) , . 

« , = e), k— ; = <0, (12) 
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where k = {k + 2)/2. It is easy to check that the zero curvature conditions 
are fulfilled: 



The solution of (O) is 



dfa 



X(v,e) 



2 



a</3 ja 



The function defined through multiple contour integrals 



(13) 



generates a solution to the (perturbed) Knizhnik-Zamolodchikov (KZ) equa- 
tion W. 



K = T,-$. 

This constitutes a generalization of the solutions of the KZ equation found 
2J 



m 



fsee also, |24 



j. [26|). It should be stressed that the perturbed 
KZ equation is understood in the sense of perturbed conformal field theory 
27[ . That is, the conformal invariance is broken once the term —2Sj/g is 



included. The above construction is an extension of the corresponding results 
for the BCS model, which is equivalent to Sierra's free field realization of the 
su(2) k WZNW model " 
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As Sierra noted for the BCS model, the on-shell 
eigenvalue of the integrals of motion of the model follows from the saddle 
point approximation in the singular limit k = —2. 



6. Scalar products and norms. 



There is an elegant result due to Slavnov |29|, [30], |3TJ] which gives a deter- 
minant representation for the scalar product of an eigenstate of the transfer 
matrix @ with an arbitrary state. In the quasiclassical limit, the leading 
terms in rj of Slavnov's scalar product gives rise to the scalar product of the 
states of this model yielding the result 



(w, M\v,M) 



n 



/3<a 



Wr 



-det M J({v a }, {wp}), 
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where the matrix elements of J are given by 

c 

v b -w b 



j ab = ^ j; 



" E 7 Tt • ( 14 ) 

^ iy a - v a )(w b - v a ) I 

Here {v a } are a solution to the equations (§), whereas {w/?} are arbitrary- 
parameters. Specializing {«;&} = {u a } for a solution to @ gives the square 
of the norm of the corresponding eigenvector @, 

(v, M\v, M) = det M K, (15) 

where the matrix elements of K are given by 

" 2 ^ 2 



Kaa E r« _*.V! E 



= 7 vr ( 16 ) 



7. Correlation functions and form factors. 

For any operator x and state \w,M) the correlation function is defined 

(w,M\ X \w, M) 



C( X ,w,M) 



(w,M\w, M) 



Of particular interest in the present case are the correlation functions for the 
intergrain current operator 



i[H, E%(!) -«i(2)] 
2f^^(6j(l)6 fc (2)-6j(2)6 fc (i; 



where W + = 6 f (l) - 6 f (2), 1¥ _ = 6(1) - 6(2). In the case when \w, M) is an 
eigenvector of the Hamiltonian it can be demonstrated that C(j,w, M) = 0. 
For general states however it is non-zero. In order to evaluate the expectation 
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value of the current for a general state requires knowing the off-diagonal form 
factors. 

We can use flTiD to obtain expressions for off-diagonal form factors for 
the spin operators. Observe that the inverse problem for the operators S~ is 
easily solved in the present case through 

S m = \im (u-e m )B(u). (17) 

A consequence is that the following off-diagonal form factors 
(w,M\S+\v,M + 1) = (v,M +l\S~\w,M) 
can be computed directly from (|14D . We find the following expression 



(w,m + i\s-\v,m) = nyW^) detM +1 T( m W,K}) ( ) 



with the matrix elements of T given by 

M+l 



i (Vb-ej)(wa-ej) 



Oi — 1 \j = 

M 

-Y.-, w r), 6<M + 1, 



^,M+l(w) 



The problem of obtaining analogous expressions for the operators W~ 
is not as straightforward and we will return to this problem at a later date. 
We note that general correlation functions and form factors can be computed 
using the leading terms in r) (quasi-classical limit) of the results derived by 
Kitanine [32] for the spin 1 Takhtajan-Babujian model. 



8. Conclusion. 



We have proposed a model describing strong Josephson tunneling in cou- 
pled small metallic grains which are described by reduced BCS Hamiltonians. 
At a particular value of the Josephson coupling energy we have demonstrated 
that the model becomes integrable, which allows for the exact calculation of 
the energy spectrum and form factors of the model. We believe that these 
results will provide insight into the nature of the Josephson effect at the 
nanoscale level where mean field approaches are not applicable. 
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